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Abstract. The breaking of the taste symmetry is studied in the temperature range between 140 MeV to 
550 MeV. In order to investigate this violation we have calculated the screening masses of the various taste 
states fitting the exponential decay of the spatial correlators. The computation has been performed using 
dynamical Nf = 2+1 gauge field configurations generated with the p4 staggered action along the Line of 
Constant Physics (LCP) defined by a pion mass of approximately 220 MeV and the kaon mass tuk 
equals 500 MeV. For temperatures below the transition an agreement with the predictions of the staggered 
chiral perturbation theory has been found and no temperature effect can be observed on the taste violation. 
Above the transition the taste splitting still shows an 0(a 2 ) behavior but with a temperature dependent 
slope. In addition to the analysis done for the pion multiplet we have performed an analogous computation 
for the light- strange and strange mesons and also looked at the scalar, vector and axial vector channels 
to understand how the multiplets split at finite temperature. Finally the temperature dependence of 
the pion decay constant is investigated to get further information about the chiral symmetry restoration. 
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1 Introduction 

One of the most famous problems related to the lattice 
discretization of quantum chromodynamics (QCD) is the 
presence of unphysical fermionic modes called doublers. 
Several methods have been proposed in the literature to 
avoid this proliferation and to reduce their number. One 
of these proposals is the well known and widely used stag- 
gered fermion formulation introduced for the first time in 
PP by Kogut and Susskind. Indeed after the staggering, 
the sixteen doublers reduce to four species that are usu- 
ally called tastes to distinguish them from the ordinary 
flavors. 

This construction has become popular in simulations of 
lattice QCD, in particular at finite temperature because 
it is computationally cheaper and preserves a remnant of 
the chiral symmetry. However, since the doubling problem 
is not completely solved, to reduce further the tastes and 
restore the correct number of degrees of freedom a rooting 
procedure has to be introduced. The validity of this tech- 
nique has long been debated in the literature 2,3,4,5,6, 
I71I51I MTU1ITT11 1 2 H 13"] and even if a definite conclusion is still 
missing, numerical evidence seems to confirm its correct- 
ness namely that the non-local terms appearing at finite 
lattice spacing become irrelevant in the continuum limit. 
At non-zero lattice spacing, due to the fact that for stag- 
gered fermions the spin and the taste degrees of freedom 
are mixed, a breaking of the continuum taste symmetry 
occurs and as a consequence only one of the sixteen pi- 



ons is a true Goldstone boson. The other fifteen states 
split according to their transformation properties under 
the symmetry group preserved by the staggered formula- 
tion [mfTSIITB] . Since these taste violations contribute to 
the cut-off dependence of physical quantities computed on 
the lattice, several ways have been introduced and applied 
in the literature to reduce them and to improve the stag- 
gered action. In fact, in the p4, aqstad, stout and HISQ 
actions this improvement was obtained by using some vari- 
ants of the gauge-link smearing [Tg l lTS PP T I I^] . 
In order to have a better comprehension of the impact of 
these cut-off effects on physical obscrvablcs it is important 
to understand how taste symmetry is violated at finite 
lattice spacing. The concrete way to do that is usually to 
look at the so called taste splitting defined by the differ- 
ence between the square of the mass of a non-Goldstone 
pion and the Goldstone one = m?^ — m| 5 E|. At zero 
temperature the situation is quite clear. All the results, 
see e.g. [T5Tri)lfTT] for the most recent investigations, con- 
firm that in the pion sector the taste splitting goes like 
a 2 where a is the lattice spacing. Remarkably this behav- 
ior can be well described with a tree-level computation in 
staggered chiral perturbation theory (SxPT) [2l?ll2"4"ll2"5] . 
Indeed in the staggered chiral Lagrangian taste-breaking 
four-fermion operators appear, physically due to the fact 
that the exchange of gluons with momenta q ~ ir/a can 

1 In the spinor-taste basis we will indicate with 7 and £ the 
spin and the taste gamma matrices respectively. 
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change the taste of the staggered quarks. 
At finite temperature the situation is not so well under- 
stood. S%PT is formulated at zero temperature and can 
not help in the understanding of taste symmetry breaking 
because in principle the structure of the chiral Lagrangian 
could be modified at finite T < T c , where T c is the tem- 
perature of approximate restoration of chiral symmetry, 
possibly leading to temperature effects. For temperatures 
above the transition we expect the breakdown of S^PT. 
The aim of this investigation is thus to go further into the 
problem and to study the dependence of taste breaking on 
lattice spacing as well as temperature. In order to extract 
this information we will calculate the screening masses of 
the different taste components of the pion multiplet at fi- 
nite temperature in 2 + 1 flavor QCD using an improved p4 
|18j staggered fcrmion action at bare quark masses which 
correspond to a Goldstone pion mass of 220 MeV and a 
kaon mass of 500 MeV at zero temperature. The temper- 
ature will range from about 140 MeV to 550 MeV and our 
study will be performed at two different sets of the lattice 
spacing corresponding to N T = 6 and 8. Our analysis will 
not be restricted to the pion sector, we have studied also 
the scalar, the axial- vector and the vector channels to ob- 
tain useful indications about taste violations also in those 
channels. Additionally we address the effect of the flavors 
on the taste breaking, analyzing some local and non-local 
pseudoscalar operators in light-strange (us) and strange- 
strange (ss) channels. 

The screening mass analysis is relevant not only for the 
understanding of the taste symmetry breaking but also 
for the study of the chiral symmetry. It is a well known 
fact that an advantage of the staggered formulation is to 
preserve a U(l) e subgroup^ of the chiral symmetry in the 
continuum. Since this subgroup relates different classes of 
mesonic operators, in the symmetric phase it predicts the 
degeneracy between certain channels. Above the transi- 
tion its restoration demands that at zero quark mass 

c(z) AV = (-iyc(z) v (i) 

where Cy and Cav are the correlators in the vector (7^) 
and axial vector (757^) channels respectively. 
Finally there is a general motivation to study screening 
masses since from them one can extract information about 
some non-perturbativc properties of the quark - gluon 
plasma. Indeed the inverse of the screening mass describes 
the distance beyond which the effect of a test hadron put 
into the medium is screened. 

The plan of the paper is the following: after an introduc- 
tion to staggered mesons (section [2]) where we briefly dis- 
cuss the symmetry properties of the staggered formulation 
both at zero and finite temperature, we present a summary 
of the predictions of staggered chiral perturbation theory 
on the taste breaking (section [3]) . Section (0| is devoted 
to the presentation of the results regarding the screening 
masses of the different mesonic states at finite tempera- 
ture before in we will summarize our results, discuss 

2 This subgroup should not be confused with the anomalous 
U(1)a of continuum QCD. 



some open problems and suggest possible directions for a 
future investigation. 

2 Staggered Mesons 

The symmetry properties of the system are fundamental 
to the analysis of the correlation functions and the deter- 
mination of the screening masses. In fact all the states can 
be classified according to their transformation properties 
under the subgroup preserved by the lattice discrctized 
theory. For these reasons in this section we review some 
basic and known facts about the classification of the stag- 
gered meson correlators starting with zero temperature 
[T4ll2"8] and then passing to the finite temperature case 

In the continuum limit the symmetry group of one stag- 
gered fermion is SU(2) x £/(4) with the first factor as- 
sociated to the rotational symmetry and the second one 
with the taste symmetry. At non zero lattice spacing this 
invariance group is broken down to the isometries group 
of the lattice. Since we are looking at zero momentum 
states with a fixed quark number, staggered mesons can 
be classified according to their transformation properties 
under a subgroup of the RF groups, the so called geomet- 
ric time slice group (GTS) that is generated by the spatial 
inversion 7 S , the taste transformations and the rota- 
tions of 7r/2 angle R kl . The staggered quark field x( x ) de- 
fined at the point x = (x, y, z, t) on a lattice with spacing 
a forms an 8 dimensional irreducible representation (8- 
irrep) of the GTS. Under the subgroup of the discretized 
version of the cubic rotation group Oh it decomposes as 
© A{ 8 F£ © Ff where a\ and Ff are, respectively, 
the trivial and the three dimensional vector representation 
with parity ±. To construct the meson states one has to 
take the tensor product of two staggered fields and just by 
counting one obtains sixty-four states that, following [14] . 
we list in Table [TJ At a first sight the meson operators can 
be divided in four classes: the local operators (n = 1..4 in 
Table[l| written as 

M°- L = 0(x)x(x) X (x), (2) 

with 0(x) being a phase factor depending on the choice of 
the channel, where the two staggered fields sit at the same 
lattice point; the one-link operators (n = 5. .10) written 
schematically as 

M l ~ L = 0(x)x(x)Ax(x) (3) 

where one has to introduce the shift operator as 

Ax(x) = l/2(x(x + ?)+ X (x-?)) (4) 

and put the quark and the anti-quark in two next-neighbor 
lattice points; and finally the two and three link operators 
(n = 11. .16 and n = 17. .20 respectively) defined as 

M 2 ~ L = e^ fc 0(x)x(x)^Zi J -x(x) (5) 

3 The rest frame group is the invariance group of the transfer 
matrix 
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M 



3-L 



(x)x(x)Z\!Z\ 2 Zi3X(x) 



(6) 



where quark and anti-quark are separated by two and 
three links. More specifically, all the states can be clas- 
sified according to the value r""" 123 of the representation 
of GTS in which they lie, where r is the dimension of the 
irrep., and a s and 0123 are, respectively, the eigenvalues 
of I s , the spatial inversion, and of X123, the parity under 
spatial taste inversion. Just for completeness, in the fourth 
column of Table Q] we have also indicated the meson de- 
scription in the original language of Kogut and Susskind: 
using the spin (g) taste basis the staggered mesons can be 
written as M = i/S (T D £g> r T ^j ip where ip has four Dirac 
and four taste components and is defined on a coarse lat- 
tice with spacing 2a. From the fifth column one can read 
off the Jr quantum numbers of the lowest correspond- 
ing continuum states at zero temperature and finally in 
the last column the particle identification of those states 
is given. 

At non-zero temperature the situation is different due to 
the fact that the temporal direction is special since it is re- 
lated to the temperature through the relation T = 1 /aN T 
where N T is the number of lattice points in that direc- 
tion. For spatial correlation functions, which we are go- 
ing to analyze, the subgroup of Oh preserved in the zero 
temperature case is now broken since certain 7r/2 rota- 
tions (for example in the x — t and y — t plane) are no 
more symmetries of the lattice discretization. As a conse- 
quence, the mesonic states have to be classified differently 
from Table [TJ The new invariance subgroup is D\ , the 
discrctized version of the cylinder symmetries with 16 ele- 
ments divided into eight one-dimensional representations 
labeled by , , B^ , and two two-dimensional ir- 
reps E . Following [22] it is easy to understand how the 
irreps of Oh reduce under D\. Basically the one dimen- 
sional irrep of Oh continues to be a 1-irrcp of D\ while 
the vector of Oh decomposes into a singlet plus the two 
dimensional representation E. One can now do the tensor 
product of two staggered quark fields and find in which 
way the mesons transform under the invariance group at 
finite T. This has already been done by S. Gupta in [29] 
and summarizing his results we have that the three dimen- 
sional vector irreps, which in Table [1] are labeled by the 
numbers n = 3. .9, 11. .14, and 19, 20, split into a two plus 
a one dimensional representation. Concerning the tensor 
irreps. of dimension 6, labeled by n = 9, 10, 15, 16, they 
split in four plus two dimensional representations. 
If the temperature is further increased and approaches the 
infinite temperature limit, compared to the finite temper- 
ature case, more degeneracies will develop as in this limit 
the theory approaches the free case. Indeed to the extent 
that dimensional reduction is a good approximation, one 
sees an approximate symmetrjjf], thus the represen- 
tations A\ , E + , A{ become degenerate with A~^ , E~ , A\ 
respectively and as a consequence the correlators fall into 
only very few multiplets. 

In this work we study the local and non-local mesons built 
from ud, us and ss flavor combinations. In flavor singlet 
channels like e.g. ss only the quark-line connected part 



has been computed and the computationally demanding 
disconnected part has been neglected. The connected part 
of the staggered meson screening correlators, projected to 
zero transverse momentum, = (Px,Py), and to zero 
(boson) Matsubara frequency w n , are obtained as 

C ( 2 ) = E E ^ (( A '4~+v,otMx)) (M-^-w))*) , 



■x,y,r v.wGl 



(7) 

where M ^ is the full staggered quark propagator from 
to x, 0(x) is a phase depending on the choice of the 
channel and £ = {0} for the local correlators, £ = {±i} 
for the one-link correlators and is equal to {±i ± j} and 
{±i±y±f} respectively for the two-link and the three 
link correlation functional. Thus while for the computa- 
tion of the screening mass of the local meson only one in- 
version of the fermion matrix is needed, for the one-link, 
two-link and three-link operators we need three, five and 
nine inversions respectively that make the computation 
more demanding compared to the local operators. Since 
a staggered fermion meson correlator, in general, contains 
two different mesons with opposite parity |14j one has to 
parametrize it as 



C(z) = Amo cosh 



M_ z- 



(-l) z A cosh 



M, z- 



(8) 



According to our phase convention, for the local and two- 
link operators M_ (-M+) corresponds to the screening 
mass of the lightest negative (positive) parity state and 
comes as the non-oscillating, NO, (oscillating, O,) part 
of the screening correlator and vicevcrsa for the one and 
three-link operator^. 



3 Staggered Chiral Perturbation Theory 

In this section we review some elements of staggered chi- 
ral perturbation theory ( SxPT) and more in detail its 
prediction regarding taste symmetry violation. The chiral 
lagrangian of the S%PT has been introduced in [23] for a 
single staggered flavor and generalized to multiple flavors 
in [2"4"ll2"6"] . It describes the low energy dynamics of pseudo- 
Goldstonc bosons near the continuum and chiral limit and 
provides a systematic method to discuss the consequences 
of chiral symmetry breaking in the framework of the stag- 
gered formulation of lattice QCD. 

Starting with an SU(3) gauge theory with n flavors of 
staggered fermions, in the combined chiral-continuum limit 
the theory possesses an 5J7(4n)i x SU(An)fj chiral sym- 
metry which is spontaneously broken to SU(An)y. The 
16n 2 -l mesonic fields that arise from the breaking can be 
collected into a unitary matrix 



E = exp 



(9) 



4 The symmetry group of the dimensionally reduced theory 
Cy ~ D4 is the group of the isometries of a square. 



5 i and ] € {£, j/,r} 

6 Note that both amplitudes are positive [31], Ano,Aq > 
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Operator 




r D x r F 
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Particle 


1 
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1++ 


1® 1 
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.1 


TV 


2 




1+- 


74 <8> £4 


or 










75 <S> £5 


or 


TV 


3 




3+" 


7i75 ® £i^5 


i ++ 


a 1 








7i74 ® £i^4 




p 


4 




3++ 


7j7fc ® £j£ib 


ir 


bi 








7i <8> £i 




p 


5 


XViAiX 


3-+ 


7< ® 1 


!?" 










7i7fc (8 £4^5 


ir 


bi 


6 


ViCiXViAiY 


3"" 


7*74 55 £4 


I4" 


p 








7i75 ® £5 


ir 


a 1 


7 




3~~ 


75 ® £i£s 


or 


7T 








74 ® C»C4 


or 




8 


/ 1 TV* /V 1 /V- 


3"+ 


7475 ® £i£fe 


0a + 


7T 








1 ® fi 


or 


do 


9 


^ieCiXCi^i'X 
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7fe74 ® £i^5 


I4" 


f) 








7fc75 <8> £i^4 


A 


"I 


10 


V4, CiVi ed XCiAjY 


6"+ 


7fc ® fiffc 


A 


p 








7i7i 55 fi 


ir 


bi 


11 


£nkXViAi(r)iAj'x) 


3++ 


7i7j 55 1 


4" 
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p 


12 


£ijkViC4XVi^i{Vj Aj'X) 
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u + 
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7fc74 55 ?5 
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P 


13 
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or 
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7475 55 £fc 
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14 


eijfe??4C4XCi /i »(0A'X) 


3+- 


74 55 CfeCs 


or 


- 








75 ® Cfe^4 
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7T 


15 


rjfcCfcX^iACCiAx) 


6++ 


7i7fe 8) 


ir 


bi 








7j 55 & 




/' 


16 


7?4C4??fcCfcX^^ l (CjAx) 


6+- 


7j7s 55 


1 A 


«i 








7j74 55 &$4 






17 


X??l^l('72^2(?73^i3X)) 


1-+ 


7475 ® 1 


or 


'/ 








1 55 


or 




18 


?74C4X ) 7l- 4 l( ? ?2^2(7?3^3X)) 


1"" 


75 55 £4 


0a + 


7T 








74 55 75 


or 




19 


»7iCCiX»7l^l('72^2(»?3^3X)) 
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7i74 55 Ci?5 
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7i75 55 £i?4 
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Table 1. Irreducible representations and their operators for mesonic states. Here Aif(x,t) = | [/(a; + i, t) + /(# — while 
»j M (ic) = C M (a;) = (-l)^+i+-^ and c(x) = 



with 



/ C/ 7T+ A'+ ■ ■ ■ \ 
7T- D A° • • • » 

A"" A S 



V 



/ 



where U = E„=i tt+ = ^ali 7T+T ai etc. are 4 x 4 

submatrices expressed in the Hermitian basifl 



7 Cm are the Dirac gamma matrices and £1 = I4 X 4 



(10) 



and /„• is the pion decay constant. Using the usual power- 
counting scheme to derive the staggered chiral lagrangian 



P 2 /Al 



l QCD ~ m I^QCD ~ a A QCD 

one obtains to lowest order 



(11) 



C 



S X PT 



f 2 11 f 2 

^Tr(d„Sd^) - ^-Tr{M^S + ME*)- 
8 4 
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-ml (U 2 + D 2 + S 2 + ---)+a 2 V(U) 



(12) 



where 



M= 



/ m u I • • • \ 
m d I •• 
m s I ■■ 



\ 



J 



is the quark mass matrix, mo models the effect of the 
anomaly and V(S) is the taste-breaking potential. This 
potential arises from the four-fermion operators in the 
quark effective action and is a linear combination of oper- 
ators 



V{S) = Y,G i O i 



(13) 



with Ci six unknown low energy effective constants (LECs) . 
More in detail the operators that enter in the potential 
read as 



(h 



Oi = Tr ( -s^Ssr 1 ") 
^[Tx(S ll E)lr{5 )t E) + h.c.] 
° 2A = I E [ Tr ( S ^E) Tr {S^E) + h.c] 
3 = i^[Tr(H M 17H M r) + ft. C .] 
4 = ^[7r(S l > 5 ZS fl5 E) + h.c] 
° 5V = ^ E [ Tr Tr ( H^t) + 

° 5A = ^ E [ Tr (S^Tr (H MS rt) + 

06 = ^E rr {s^es^ r f ) 



/i.e. 



/i.e. 



(14) 
(15) 

(16) 

(17) 

(18) 

(19) 

(20) 

(21) 



Here the An x An matrices St are just the generalizations 
of the 4x4 taste matrices £t- If we focus only on the 
non-diagonal flavor meson^f] it is sufficient to expand the 
lagrangian to quadratic order in the mesonic fields and 
consider only the single-trace terms. ^From the expansion 
one can read off the tree-level masses of the mesons : 



m M B = K m a + m b ) + a 2 A^ E 



(22) 



where the meson M is composed of two quarks a and b 
and where 



A(&) = A PS = 



(23) 



Simulations in which disconnected contributions are not 
taken into account describe only non-diagonal flavor states 
n ± ,K ± .. 



A(^ 5 ) = A A = ^(C 1 +3C 3 +C 4 + 3C 6 ) (24) 

J 7T 



A(^) = A T = ^(2C S + 2C 4 + AC 6 ) (25) 



^ (Ci + C 3 + 30 4 + 30 6 ) 



16 

■' 7T 



4C 4 ) 



(26) 



(27) 



For the flavor neutral mesons (U, D, • • • ) the situation is 
more complicated since in the expansion of the two-trace 
operators some two-point vertices mixing the taste- vector, 
taste-axial and taste-singlet flavor-neutral states appear 
and thus the Lagrangian receives a contribution of the 
form 



C 



16a 2 , 

—S V (U^+D^ + 



f 

16a 2 
f 2 

J 7T 



+ 



SavP^+D^ + ■■■)* 



J 7T 



(28) 



For more details about the taste splitting in such channels 
we refer to [2~4] . 

One interesting thing to note is that the symmetry group 
of the lattice theory at 0(a 2 ) is enlarged in the meson 
sector since the potential V(S) is rotationally and 50(4) 
taste invariant. The consequence of this enhancement of 
symmetry is that the sixteen taste states instead of falling 
into eight different irreps of the GTS as predicted from 
group theory split into only five degenerate multiplets ac- 
cording to the 5*0(4) symmetry (eq. l2"3"ll2~7*j) . The vanishing 
of A(£s) is due to the taste non-singlet U e (l) symmetry 
which is unbroken by the lattice regulator, making the 
75 x £5 meson a true Goldstone boson. This prediction of 
the splitting has been well confirmed in a series of lattice 
simulations [TSJUMnil^ ■ 

This analysis is strictly valid at temperature zero. Since 
our aim is to study the taste symmetry violations at finite 
T we need to understand how the predictions of S%PT 
get modified. In particular two scenarios seem possible. 
The first one is that the formal structure of the chiral La- 
grangian remains the same as in (|T2l but the coefficients 
Ci and /,r acquire a temperature dependence. At least 
this hypothesis could work until T c where the chiral sym- 
metry is restored and chiral perturbation theory breaks 
down. Indeed, as we will see, in the vicinity of the transi- 
tion the value of the pion decay constant that " measures" 
the strength of the chiral symmetry breaking goes to zero 
rapidly. A second scenario can be described by a different 
lagrangian since in principle operators different from the 
T = case are admitted in the Symanzik effective theory. 
This is due to the fact that, as explained in the previ- 
ous section, the invariance group of the zero temperature 
theory is broken down to a subgroup. 
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775 
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0.00253 
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0.00259 


240 
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3.530 
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590 


3.540 


0.00240 


259 
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3.570 


0.00212 


281 


525 
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738 


o.ooo 








91 Q 


ouu 


3.630 


0.00170 


326 


389 






3.690 


0.00150 


365 


448 






3.760 


0.00130 


424 


465 


318 


588 


3.820 


0.00125 






361 


488 


3.920 


0.00110 






410 


481 


3.920 


0.00092 


532 


480 






4.000 


0.00092 






475 


493 


4.080 


0.00081 






549 


392 



Table 2. Coupling constants /3 and light quark masses mi, 
temperatures T in MeV and the number of configurations of 
the given sizes on which screening masses were computed. The 
strange quarks have always been choice as m 3 = 10m;. 



4 Results 

In order to perform this analysis we have used dynami- 
cal Nf =2+1 gauge field configurations generated with the 
RHMC algorithm 02 by the RBC-Bielefeld [HMS] and 
the HotQCD [M] collaborations using the p4 staggered 
action. The configurations were generated along the Line 
of Constant Physics (LCP) obtained by tuning the bare 
quark masses such that at zero temperature the (Gold- 
stone) pion mass is approximately 220 MeV and the 
kaon mass ttik equals 500 MeV. The bare strange quark 
mass has always been 10 times larger than the light one, 
m s = 10m;. In table ©, for both lattices 24 3 x 6 and 
32 3 x 8 used in the simulations, the values of the cou- 
pling constant and of the bare light quark mass at which 
the configurations were generated are reported. The num- 
ber of configurations analyzed was about 500-600 at each 
temperature, separated by 10 time units. For more details 
about the simulations and scale settings we refer to |33j . 
In the following we report on our results for the screen- 
ing masses of the local and non-local staggered mesons 
extracted from the lattice calculation of the spatial cor- 
relators in some of the twenty mesonic channels listed in 
Table ([T]). More in detail we focus on the pion channels 
labeled with n = 1,2, 7, 8, 13, 14, 17, 18 since we are pri- 
marily interested in the taste splittings in this multiplet. 
However our analysis is not limited to that but we also 
inspected some scalar, vector and axial- vector channels to 
extract some other useful information from them, for ex- 
ample about the restoration of the chiral symmetry and 
the finite temperature splittings of those multiplcts. 



Before starting the presentation of the numerical results 
we recall that a staggered correlator contains two differ- 
ent mesons with opposite parity (see eq. . This is valid 
in general except for the channels 2, 7, 14, 18 where only 
the pion contributes. From these channels we are able to 
extract the screening masses above and below the transi- 
tion fitting as usual the correlators with the exponential 
functions and with errors determined by jack-knifing the 
fit. In the other channels the pions have to be extracted 
from the oscillating (n = 1, 13) or the non-oscillating part 
(n = 8, 17) but also the scalars contribute to the corre- 
lators. If one performs the analysis at zero temperature 
no problem occurs in the determination of the screening 
masses of all pion taste components. Unfortunately the 
situation is quite different at finite temperature since in 
the channels n = 1,8,13,17 the amplitudes of the pion 
states die out very fast with rising temperature and as 
a consequence we are not able to extract their screening 
masses. Indeed at temperatures above the transition the 
pion contributions to the correlators disappear completely 
and the channels are dominated by the scalars. The lack 
of these screening mass data makes the comparison with 
the S%PT incomplete. Differently from the T = case 
in which one can check an SO(A) restoration of the taste 
symmetry predicted by the SxPT, we are not able to verify 
if an analogous enhancement happens also at finite tem- 
perature. 

For the axial-vector and vector channels the situation is 
not so clear, we have seen in all analyzed channels the 
presence of both parity states but, likewise to what hap- 
pens in the local sector |35j . using the point sources we 
were able to extract only few values of the screening mass 
for very high temperatures (above 2 T c ) since the correla- 
tors are very noisy. 

In addition to the investigation done for mesons made out 
off two degenerate u and d quarks we have also performed 
some calculations for meson with flavor content us and 
ss to compare the result with the light sector and extract 
some additional checks regarding the S%PT that predicts 
the flavor independence of the taste violation. 
Let us start with figure []] where we present the data on the 
screening mass of the various taste components of the pion 
multiplet labeled by their taste matrix [£]. In particulari, 
in the pictures l(a)| and 1 (c) we plot the screening masses 
in units of rg as function of the temperature at N T = 8 and 
N T = 6 respectively. For both lattices we can observe that 
below the transition the Goldstone boson screening mass 
([£] = £5) remains approximatively constant while for the 
non local operators a slight decrease of their masses occurs 
below T c . Above the transition there is a rapid linear rise 
of the values of the masses for all the states with increas- 



ing temperature. From the figures 1(b) and 1(d) where 
the screening masses in unit of T as function of the tem- 
perature are plotted one can note a strong decrease for 
the non-local operator masses below the transition and an 
increase above that point with an approach from below 
to the free continuum result given by 27rT (indicated by 
the continuum line) where the taste symmetry is restored 
completely. 
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Fig. 1. Temperature dependence of the pseudoscalar screening mass (75 ® [£]) in units of ro and of the temperature at N T — 8 
(a, b) and N T = 6 (c, d) for different components of the multiplet identified by the taste matrix [£] . 



The scalar multiplet is shown in fig. [5] where we plot the 
temperature dependence of the screening masses of the 
various taste components normalized by ro ( 2(a)| 2(c) I 
and by the temperature ( 2(b) [[2(d)) . Below the transition 
it is difficult to obtain quantitative results since the corre- 
lators turn out to be very noisy in this region, especially 
for the non-local observables. At temperature above 1.2 T c 
all scalar correlators are degenerate with the pseudoscalar 
ones in corresponding taste channels. 
In figure [3] we show the flavor dependence of the screen- 
ing mass for the local and the 1-link pseudoscalar channel. 
More precisely, we have calculated and plotted the screen- 
ing masses in the light (ud), light-strange (us) and strange 
(ss) sectors. While above 1.5 T c for both, local and 1-link 
operator no difference between the three sectors can be 
observed since the quark mass term becomes irrelevant, 
below T c Pion, Kaon and <fi split as expected. Further- 
more, in the region below the transition we cannot identify 
any substantial temperature effect on the local operator's 
masses while for the one-link operators a slight decrease 
can be observed as the temperature is increased until T c . 
Note however that this data has been obtained at fixed 



In the following we will discuss the taste splittings which 
arc defined as the differences of squared pseudoscalar masses 
of different taste, see Eci.([22"j). 



A t = 



(29) 



and are traditionally but somewhat arbitrarily normalized 
to (200 MeV) 2 . In order to understand the effect of the 
temperature on the taste symmetry violation in the group 
of figures (Q| we have plotted the taste splittings as func- 
tion of T at fixed /3 values. The results indicate that while 
in the confined phase there is essentially no effect of the 
temperature on the taste splitting, as one can figure out 
from figure 4(a)[ above the transition a temperature effect 
can be clearly read from 4(b)[ 4(c) 4(d) and 4(e) This un- 
expected behavior above T c deserves some comments that 
we postpone to the discussion section. Here we underline 
the fact and note that this effect seems independent of the 
lattice spacing since the slope of the increase stays approx- 
imatively constant for the configurations at different beta 
values. 

In the group of figures [5] we study the lattice spacing de- 
pendence of the taste violation at fixed temperature, plot- 
ting the value of the taste splitting as function of a 2 . For 
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Fig. 2. Temperature dependence of the scalar screening mass (1 fg> [£]) in unit of ro and of the temperature at N T = 8 (a,b) 
and N T = 6 (c,d) for different components of the multiplet identified by the taste matrix [£]. 
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Fig. 3. (a) Flavor dependence of the screening mass for the local pseudoscalar operator (75 ® £5) in units of ro at N T = 6. (b) 
Flavor dependence of the screening mass for the 1-link pseudoscalar operator (75 ® ^£5) in units of ro at N T = 6. 
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Fig. 4. Effect of the temperature on the taste violation at fixed /3 values : (a) j3 = 3.43, (b) /? = 3.50, (c) /3 = 3.57, (d) j3 = 3.76, 
(e) /3 = 3.92 



all analyzed temperatures the taste splittings arc compat- 
ible with an a 2 behavior. However as we can easily see 
from 6(a) and |6(b)] while below the transition the slope 
is more or less constant above T c its values start to in- 
crease with the temperature. This is an indication that in 
the low temperature regime the predictions of the S%PT 
on the taste splitting seem to be valid and independent 
from T; even if above the transition these predictions are 
no more reliable an a 2 behavior is still observed but the 
temperature starts to affect the taste violation making the 
quantity A^] temperature-dependent. 
In figure 7(a) we have plotted the taste splitting for dif- 
ferent taste states in units of 200 MeV as function of the 
temperature. The dependence of the taste violation can- 
not be extracted trivially looking at this pictures since 
in the plot we are comparing different finite temperature 
configurations with the same N T but thus with different 
values of t he lat tice spaci ng. For example comparing this 
figure with [6(a) and 6(b) one can note that the stro ng de- 
crea sing o f the taste splitting observed below T c in (7(a) I 
and ( |7(b)[ ) is just a pure discretization effect. Above T c the 
behavior of the splittings changes drastically, first start- 
ing to rise rapidly and at 1.5 T c there is a flattening of 
the slope but again this is a combined lattice spacing and 
temperature effect. 

The flavor dependence on the taste violation is analyzed in 
the figurc |7(b)| where the quantity A ia = -ro 2 5075 



for the ud, us and ss channel is plotted. In the range of 
the temperature analyzed the screening masses in the dif- 
ferent sectors are in agreement within the error bars. But 
while below the transition where the quark mass is rele- 
vant this fact is a not-trivial prediction of the staggered 
chiral perturbation theory, above 1.5 T c the plot tells us 
simply that as noted in fig. 3(a) and |3(b)| thc effect of the 
quark mass is negligible. 

In order to study the restoration of the chiral symmetry 
and the breakdown of the chiral perturbation theory it's 
interesting to study the temperature dependence of the 
pion decay constant / T that is defined on the lattice [55] 

as 



1 mi sinh (m B 



1/2 



4 sinh (m„/2)' 



PS 



(30) 



where and Aps0 are respectively the mass and the 
amplitude extracted from the Goldstone pion correlator. 
This relation is nothing else than the staggered lattice ver- 
sion |36j of the continuum definition through the coupling 
of the axial current to the pion 



V^f^m^ = -(0|u747 5 d|7r + ) 



(31) 



The amplitude Aps is related to A N q(o) (equation [SJ) 



trough the relation A 



JVO(O) 



= 2A 



pse 
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Fig. 5. Effect of the lattice spacing on the taste violation at fixed temperature values 
(c) T = 226 MeV, (d) T = 361 MeV. 
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From figure 8(a) one can essentially note the sharp de- 
crease of the value of /„■ in the proximity of the transi- 
tion and a near vanishing after that point. This signals 
clearly that the chiral symmetry is restored and the chiral 
perturbation theory is no more applicable. In the chiral 
symmetry broken phase even if only few data are avail- 
able we can observe a temperature dependence of / W (T). 



Another thing to point out is the difference between the 
N T = 8 and the N T = 6 results. This is probably due to 
the fact that in the first case the discretization errors are 
reduced compared to the second lattice and indeed the 
values are closer to the continuum value of f v (Q) — 93 
MeV. It would be interesting to get some more data in 
the confined phase and to perform an extrapolation to 
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Fig. 7. (a) Effect of the temperature on the taste violation in the pion channel at 7V T = 6. (b) Flavor dependence of the taste 
violation for the light-light, light-strange and strange-strange pseudoscalar channels. The quantity — w^ 6 ^ x — m^ 5 g, 75 is 
plotted for different flavor content at N T = 6 




Fig. 8. (a) Temperature dependence of the pion decay constant / w in MeV calculated at N T =6. (b) Temperature dependence 
for different taste state in the vector channel identified by the spin x taste matrix y x ® [£] calculated at N T =6. 



the physical mass point and to the continuum similar to 
that done at zero temperature in (38j . 
In the last picture we plot the temperature dependence of 
the screening masses for some taste states in the vector 
channel. Unfortunately since at low temperatures these 
correlators are very noisy it is in general difficult to ex- 
tract the vect or screening masses for T < 2T C . Looking 
at figure 8(b) we can essentially note two facts: the first 



is that at temperatures between 2 and 3 T c the values of 
the screening masses in the vector channel are more or less 
constant and have already reached a value of 2ttT which is 
the infinite temperature value in the thermodynamic limit. 
The second point is that taste violations do not seem to 
be present in these channels since the difference between 
the various tastes are within the error bars. We have also 
looked at other vector and axial-vector channels finding 
analogous results. 



5 Summary 

In this paper we have investigated the violation of taste 
symmetry at finite temperature in 2 + 1 flavor QCD uti- 
lizing the gauge configurations generated by the RBC- 
Bielcfcld [3"0ll3"5] and the HotQCD [33] collaborations us- 
ing the improved p4 fermion action. Two different lattice 
sizes, 24 3 x 6 and 32 3 x 8, have been used to understand 
the role of the temperature in the taste breaking phenom- 
ena. We have calculated the screening masses by fitting 
the spatial meson correlators and even if we have been 
primarily concentrated on the pion channels, some other 
mesonic states built from strange quarks and some scalar, 
axial-vector and vector operators have been considered. 
The analysis of the taste symmetry breaking is important 
since this violation is related to the cut-off dependence of 
the physical quantities. At temperatures below the tran- 
sition our aim was to understand whether the predictions 
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of the staggered chiral perturbation theory on the taste 
violation are still valid, with possibly temperature depen- 
dent effective low energy constants, or whether eventually 
some new terms in the chiral lagrangian would have to 
be considered. As pointed out in the previous section for 
T < T c we have found indications that the taste splitting 
in the pion channel has an a 2 behavior as predicted by 
S^PT and no effects on the temperature can be observed 
in this region. It could be interesting to compare the re- 
sults with zero temperature taste splittings with the same 
action. 

The apparent temperature independence of the taste split- 
tings in combination with a noticeably temperature de- 
pendent pion decay constant (fig. 8(a) I implies within 
staggered chiral perturbation theory that the coefficients 
Ci of the lagrangian (|12[) acquire a temperature depen- 
dence. However, as not all taste states in the pion multi- 
plet are accessible at high temperature due to their dying 
out when the transition region is approached, the system 
of equations for the C, cannot be solved and it is not possi- 
ble to determine the temperature behavior of these LECs. 
Above the transition the situation changes drastically. As 
also indicated by the vanishing of the pion decay con- 
stant f w , chiral symmetry is restored and the predictions 
of chiral perturbation theory are no longer valid. Staying 
nevertheless with the chirally inspired definition of taste 
splittings, eq. (|22j) . this quantity measures taste splittings 
that behave like a 2 - as is expected - but with a slope which 
increases with temperature, i.e. the splittings appear to be 
temperature dependent. However, because chiral symme- 
try is restored in this region the pion ceases to be a Gold- 
stone particle and acquires a non- vanishing mass. For suf- 
ficiently high temperatures the screening masses will be- 
come proportional to the temperature. As a consequence, 
the splittings as defined in ([22]) should then behave ac- 
cording to 



regime are available. The situation may improve with the 
use of wall sources succeeding to analyze configurations at 
lower temperatures. Quite surprisingly, however, we found 
that pion triplets n = 7,8, 13, 14 do not split in this way. 
Indeed the difference between the "transversal" and the 
"longitudinal" screening masses is compatible within the 
error bars. This is a clear indication that at finite T at 
least in the pion channel, contrary to group theory ar- 
guments, an enhancement of the 5*0(2) x Z2 symmetry 
occurs and a restoration of an SO(3) symmetry can be 
observed in these channels. As already stressed previously 
we don't have access to all the taste masses, so that we 
cannot check for a complete S'0(4) restoration like in the 
zero temperature case. 

Finally we have evaluated the flavor effect on the taste 
splitting. For temperatures above 1.5 T c the splittings are 
flavor independent for the simple reason that in this region 
the effect of the quark mass is negligible. For temperatures 
below the transition the fact that no flavor effect can be 
observed is non-trivial and supports the applicability of 
staggered chiral perturbation theory, which predicts that 
the pion splittings are independent of the flavor of the 
quark constituents, in this temperature regime. 
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z%] ~ m £5 (m K] - m Cs ) ~ T(m K] - m f5 ) > a 2 T, (32) Appendix 



i.e. rise with temperature. The temperature interval in- 
vestigated, T c < T < 2T C , represents the transition re- 
gion where the pion changes from Goldstone nature to 
approaching free theory behavior such that the temper- 
ature dependence is more complicated and a linear rise 
not yet clearly identifiable. It would be interesting to bet- 
ter understand the behavior in the chirally restored phase 
and to compare for instance with perturbative predictions 
which could guide the numerical analysis and should be 
reliable at high temperature. 

Another point that we have analyzed carefully is the fi- 
nite temperature splitting of the T — multiplets [25] . 
According to group theory e.g. the triplets of the T = 
GTS split into a two plus a one dimensional irrep of the 
finite temperature GTS. As was already noted in the liter- 
ature [35] for some local axial- vector and vector channels a 
splitting between "longitudinal" and "transverse" screen- 
ing masses has been observed. Indeed, also in the other, 
non-local triplet channels (n = 5,6,11,12,19,20) these 
splittings can be seen even if the analysis is performed with 
point source operators where only few data in the high T 



In the tables[3jand|4]we have summarized the pion screen- 
ing masses for different taste states. The screening masses 
are given in lattice units and the data can be easily con- 
verted to either temperature or vacuum (ro) units by means 
of the included zero temperature results for rg/a from [33] . 
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